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We report a numerical study on Aharonov-Bohm (AB) effect and giant magnetoresistance in 
rectangular rings made of graphene nanoribbons (GNRs). We show that in low energy regime where 
only the first subband of contact GNRs contributes to the transport, the transmission probability 
can be strongly modulated, i.e., almost fully suppressed, when tuning a perpendicular magnetic 
field. On this basis, strong AB oscillations with giant negative magnetoresistance can be achieved 
at room temperature. The magnetoresistance reaches up to a few ten thousand % in perfect GNR 
rings and about a few hundred % with edge disordered GNRs. The design rules to observe such 
strong effects are also discussed. Our study hence provides guideline for further investigations of 
the AB interference and puts forward the high magnetoresistance applications of graphene. 
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Graphene and its nanostructures have recently at- 
tracted a great amount of attention for both fundamental 
researches and device applications [THS]. This is partic- 
ularly due to its unusual electronic properties such as 
the linear dispersion and the chirality of carriers mak- 
ing graphene definitely different from conventional solid- 
state materials. These properties lead to many unusual 
transport phenomena in graphene structures such as fi- 
nite minimal conductivity, Klein tunneling, or unconven- 
tional quantum Hall effect (e.g., see the review [Tj). Addi- 
tionally, graphene also possesses outstanding properties 
such as high carrier mobility [4^ and small spin-orbit cou- 
pling which make it very promising for applications in 
electronics and for use in ballistic spin transport devices. 
Various studies in this direction have hence been carried 
out (e.g., see the reviews [11 13]). 

The Aharonov-Bohm (AB) oscillations [6^ in meso- 
scopic rings are a phenomenon of particular interest and 
an elegant way to study phase coherent transport. In the 
presence of a perpendicular magnetic-field 5, the phase 
coherent trajectories of charge carriers encircling the ring 
are characterized by the phase difference = 27rBS/(j)o^ 
where = h/e and S is the area of the ring. There- 
fore, the transmission probability through the ring ex- 
hibits oscillations when tuning the magnetic field with 
period AB = (j)o/S. The AB effect was originally ob- 
served in metal rings [7j , and later in semiconductor het- 
erostructures [8] , carbon nanotubes [9l |10) , and topolog- 
ical insulators [11]. The AB oscillations have been also 
explored in mesoscopic graphene rings (see the recent re- 
view [12 ). Experimentally, clear h/e—AB oscillations 
have been observed in monolayer graphene rings [13HT6 1, 
graphene films with antidot arrays [17 , and thin graphite 
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FIG. 1: Schematic of the graphene nanoribbon rings consid- 
ered in this work. Qr, Qc, and Qh characterize the width of 
the ring, of the contact graphene nanoribbons and of the hole, 
respectively. Nh defines the length of the hole and Ns stands 
for the size of side nanoribbons along the transport direction. 



crystals with columnar defects [T8^. On the theoretical 
side, many interesting effects have been investigated and 
discussed, including the valley degree of freedom typical 
of graphene, the influence of particular device geometries 
and edge symmetries, a resonant behavior with transistor 
applications, or the interplay between the AB effect and 
Klein tunneling \T9*-^26* . However, in almost all struc- 
tures studied previously, the phase coherence was not as 
strong as expected and hence the amplitude of AB oscil- 
lations and magnetoresistance (MR) was relatively small 
even at low temperature. This will be discussed in more 
detail in this paper, on the basis of our investigations. 

It is well known that achieving a high magnetoresis- 
tance is especially crucial for applications such as high- 
density data storage or magnetic sensors and actuators 
[27] , Hence, the investigation of this effect in graphene 
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FIG. 2: (a,c) bandstructure of the contact GNR and (b,d) cor- 
responding transmission probability of armchair GNR rings at 
B = and 8 T. Panels (a,b) and panels (c,d) are for the rings 
schematized in fig. 1(a) and fig. 1(b), respectively. The wave 
vector is given in units of 7r/ao, with ao — 3ac. Parameters: 
= 23, Qh = 13, Nh = 120 and Ns = 11 in (c,d). 



nanostructures with either ferromagnetic (e.g., see the 
review [2 ) or non-magnetic contacts [28H37] has re- 
cently been an emerging research topic. For instance, 
it has been experimentally reported [28] and theoreti- 
cally demonstrated [29] that a high MR of ~ 50 % can 
be obtained at room temperature in graphene nanorib- 
bon (GNR) devices thanks to the reduction of bandgap 
induced by the presence of a magnetic field. Similarly, 
large MR has been observed in p-i-n GNR heterostruc- 
tures [30]. In ref. [31 , a MR close to 85 % at room 
temperature has been achieved thanks to the orthogonal- 
ity of the wavefunctions in metallic and semiconducting 
GNR sections. In ref. [32 , a large (- 50 %) MR has been 
experimentally shown in multilayered epitaxial graphene. 
Additionally, the low-temperature magnetotransport has 
been also studied in various works [33H37j . 

In this article, we investigate the magnetotransport 
in the rectangular GNR rings schematized in Fig.l 
and predict strong AB oscillations with a huge room- 
temperature MR. These strong effects are observed in the 
low energy regime corresponding to the first subband of 
the contact GNRs. Our calculations show that a negative 
MR exceeding 10^ % in perfect GNR rings and 10^ % in 
edge disordered ones can be achieved. We also reach the 
conclusion that it is hard to observe such strong effects 
in the rings previously studied in the literature because 
of the multisubband contribution of contact GNRs to the 
transport and/or of their inhomogeneous geometries. 

We use the nearest neighbor 7r-orbital tight binding 
model [TJ [26l [38] to compute the electronic transport 
in GNR rings under a uniform perpendicular magnetic- 
field (5-field). In the presence of the 5— field, the tight 
binding Hamiltonian H is modified within the Peierls 
phase approximation [35l [39] . The hopping integral be- 
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FIG. 3: (a,b) conductance and (c,d) corresponding magne- 
toresistance in armchair GNR rings as a function of 5— field 
for different Fermi energies. Panels (a,c) and panels (b,d) are 
for the rings studied in fig. 2(b) and fig. 2(d), respectively. 



tween nearest-neighbor atoms is hence given by tnm = 
— To exp(z^nm), where tq ^ 2.7 eV [26^ and (j)nm = 
^ f^;^A{r)dr. The vector potential A(ry = {-By, 0,0) 

is related to the magnetic field B = (0, 0, 5) by Vx A = 
B. The charge transport through the ring is computed 
using an adaptive recursive Green's function method, ca- 
pable of treating systems of arbitrary shape [40]. The 
linear conductance and the current are calculated using 
the Landauer formula 
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where fL{R){s) = [l ^ exp {{e - EpLiR)) /hT)]~^ is 
the left (right) Fermi distribution function with the Fermi 
level EpLt^ji^ and Go = 2e^//i is the quantum conduc- 
tance. The transmission probability is computed as 
T{e,B) = Tr [FLG^ri^G^t] from the device retarded 

Green's function = [E ^ iO^ - H - - ^rV^ , 
^L{R) = i (^^L - ^l)^ self energy T^l{r) defin- 

ing the left (right) contact-to-device coupling. Fi- 
nally, the magnetoresistance is defined as MR = 
[I {B) - 1 (0)] /I {B) under a finite bias or MR = 
[g (B) - g (0)] /g (B) at zero bias. 

Let us first investigate the properties of AB interfer- 
ences in the considered rings. The ring geometry is char- 
acterized by the set of parameters of Fig. 1. The width 
of the GNRs (Qc, Qr, Qh) is given in units of ac a/3/2 
and 3ac/2 while their length {N^^ Ng) is given in units 
of 3ac and acV^ in armchair and zigzag GNR rings, re- 
spectively, with ttc = 0.142 nm. In Fig. 2, we display 
the bandstructure of contact GNRs (left panels) and the 
transmission probability (right panels) of two different 
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FIG. 4: (a,b) conductance and (c,d) corresponding magne- 
toresistance in zigzag GNR rings as a function of 5— field for 
different Fermi energies. Panels (a,c) and panels (b,d) are for 
the rings having the same geometry as in fig. 1(a) and fig. 
1(b), respectively. Parameters: Qr — 22, Qu = 14, Nu = 120, 
Qc — 58 in (a,c) and 26 in (b,d), and A^s = 11 in (b,d). 
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FIG. 5: / — y characteristics (see the left axis) at B = and 8 
T and corresponding magnetoresistance (see the right axis) of 
different GNR rings. Panels (a,b,c,d) correspond to the rings 
studied in fig. 3(a), fig. 3(b), fig. 4(a), fig. 4(b), respectively. 
The Fermi energy is Ep — 0.1 eV in (a,c) and Ep — 0.2 eV 
in (b,d). 



armchair GNR rings: Figs. 2(a,b) for the ring shown in 
Fig. 1(a), and Figs. 2 (c,d) for the ring of Fig. 1(b). Both 
the contact and ring GNRs are metallic with a negligi- 
ble bandgap. The presence of a 5-field does not affect 
significantly the bandstructure of the contact GNRs be- 
cause they are not large enough. However, an interesting 
phenomenon is found: due to the AB interference (shown 
below), the transmission probability can be strongly sup- 
pressed in the energy regime corresponding to the first 
subband of the contact GNRs, while the influence of the 
5— field is weak at higher energies. We suggest that these 
features can be understood as follows. At low energy, the 
contacts inject a pure state of incoming wave into the 
ring and the AB interference can be perfectly achieved. 
At high energies, i.e. when several subbands can carry 
electrons, the incoming wave is no longer a pure state 
and hence the AB interference can not take place prop- 
erly. We find that these features can be reproduced in all 
rings with different parameters, regardless of the metal- 
lic/semiconducting or armchair/zigzag character of the 
GNRs (see below). In the general case, the energy regime 
where a strong AB interference takes place is determined 
by G EsAB = [^ei,^e2], in which Ee2 is the second 
band-egde of the contact GNRs and E^i is the lowest of 
the first band-edges of the contact and ring GNRs. This 
is a key-point, which motivates us to investigate the AB 
interference and the possibility to obtain high magne- 
toresistance in the proposed rings. Regarding the rings 
schematized in Fig. 1(b), we focus here on the cases of 
Qc > Qh (at variance with the studies in [24l |25]) to 
observe a strong MR effect, as discussed later. 

To clarify how strong the AB effect can be, we plot 
in Fig. 3 the conductance and the corresponding MR 
as a function of 5— field for different Fermi energies in 



the two rings studied above. Note that in what follows, 
all transport quantities are calculated at room tempera- 
ture. It is shown that (i) the conductance exhibits clear 
AB oscillations (see Figs. 3(a,b)), the period of which 
matches well the expression AB = (j)o/S, i.e., AB ^ 16 
T for 5 ^ 258 nm^; (ii) a giant negative MR of about a 
few thousand percents (see Figs. 3(c,d)) can be achieved. 
For completeness, we display in Fig. 4 the data obtained 
in rings made of zigzag GNRs. Similarly to the arm- 
chair cases, strong AB oscillations with giant MR are 
obtained. However, the transport at low energy takes 
place in the zigzag rings mainly via the edge localized 
states in the GNR arms, which weakens the confinement 
effects [24l|25]. Hence, the transmission probability (and 
conductance peaks) in the phase coherent cases is higher 
than that in the armchair rings. This leads to AB oscil- 
lations of large amplitude (see in Figs. 4(a,b)), so that 
an extremely strong MR of up to a few ten thousand per- 
cents (see Fig. 4(d)) can be achieved for the ring of Fig. 
1(b) with a large Esab- A similar giant modulation of 
the conductance, which was explained by the presence of 
field induced energy gap, has been also explored experi- 
mentally in ballistic carbon nanotubes [10 . 

Next, we explore the I — V characteristics of the con- 
sidered rings. In Fig. 5, we display the I — V curves of 
the four rings studied in Figs. 3 and 4. Interestingly, a 
giant MR can still be obtained in the finite bias regime. 
At variance with the devices studied in [28^-'3r where 
the conduction gap at low bias is reduced, the structures 
considered here can switch from metallic to semiconduct- 
ing behavior with an enhancement of the conduction gap 
when applying a 5— field. By changing the width of the 
contact GNRs, we can tune the value of Eqab and hence 
the bias window where the conduction gap takes place. 



4 



400 

gsoo 

^200 
100 



i 1 


1 1 


A| (a) 






























20 30 
Magnetic field (T) 

FIG. 6: Magnetoresistance at zero bias as a function of 
B— field with different edge disorder configurations. Panels 
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Note that the possibility of enlarging Esab (by reducing 
the width of the contact GNRs) is an advantage of the 
ring of Fig. 1(b) compared with that of Fig. 1(a). These 
results are very promising for the design of magnetic tran- 
sistors as proposed in [28l|3T]. Moreover, a specific fea- 
ture, the appearance of low (even negative) differential 
conductance at high bias, is observed in the zigzag rings. 
This is nothing but a consequence of the parity selection 
rule, which has been shown to result in the negative dif- 
ferential conductance in zigzag GNRs with an even num- 
ber of zigzag lines [4TJ 02] . This rule does not apply if 
this number is odd and hence the 1 — V characteristics of 
the zigzag rings with an odd number of zigzag lines (not 
shown here) have a form similar to that of the armchair 
ones. 

One more important point to consider is the effects 
of edge disorder, which are known to degrade the per- 
formance of most GNR devices. In Fig. 6, we display 
the MR as a function of 5— field with different disor- 
der configurations in the two rings studied in Figs. 3(a) 
and 4(a). The edge disorder is simply generated by ran- 
domly removing the edge atoms with a probability Pd. 
Indeed, the disorder strongly affects the results, i.e., it 
is hard to completely switch off the current with the AB 
interference and hence the MR amplitude is much re- 
duced. This is due to the fact that on the one hand, the 
electronic properties of the system are strongly modified 
by the disorder and, on the other hand, complex phase 
shifts are induced by the scattering of wavefunctions by 
the defects along the ring arms. These two effects totally 
weaken the AB interference. However, it is worth noting 
that in the disordered rings studied here, a large MR of 
a few hundred percents can still be achieved. Moreover, 



with the rapid progress in carbon technology j43H46] , one 
can optimistically expect that the fabrication of GNR 
rings with atomically precise edges can be achieved soon. 

Finally, we would like to discuss the reasons why it 
has been hard to obtain the strong AB inteference in the 
structures previously studied in the literature. First of 
all, in the rings made of large contact GNRs, the AB 
interference is relatively weak because of the multisub- 
band contribution to the transport as shown above. Ad- 
ditionally, with respect to the rectangular rings studied 
here, the other systems always suffer from strong inho- 
mogeneities along the transport trajectories, as a conse- 
quence of the irregular edges of the GNRs in the circular 
rings [T3H11 [H HO], or of the mixing of different GNR 
sections in the other geometries [21] [221 1241 125J. The 
inhomogeneities along the ring arms can result in com- 
plex phase shifts and hence weaken the AB interference, 
similarly to the disorder effects discussed above. These 
two reasons can explain the small MR obtained in the 
literature, compared to the strong effect observed here. 
On this basis, we note that in spite of having a large 
EsAB^ the ring geometry of Fig. 1(b) requires a careful 
design. Actually, there exists a mixing (with significant 
fractions) of zigzag and armchair GNRs in this ring if the 
side GNRs are too long. Therefore, the AB interfence can 
not take place properly for too large Qh/Qc ratio and/or 
short Nh. In that case, only the resonant tunneling effect 
due to the ring geometry is significantly pronounced as 



reported in 



The condition Qc ^ Qh is hence 



mandatory to guarantee the strong MR effect. 

In summary, we have investigated the AB effect in rect- 
angular GNR rings using numerical simulation within a 
tight binding model. We have shown that in low energy 
regime where only the first subband of contact GNRs 
contributes to the transport, i.e. in the case of a pure- 
state incoming wave, the transmission probability can 
be almost fully suppressed due to the AB interference. 
This suggests the possibility of tuning the structure from 
metallic to semiconducting state. Very strong AB os- 
cillations with huge magnetoresistance can be achieved 
at room temperature (up to a few ten thousand % in 
perfect GNR rings and about a few hundred % in edge 
disordered GNR rings). The difference between the re- 
sults obtained here and in the literature has been also 
discussed. The study hence suggests an efficient way 
to investigate the AB interference in graphene nanorings 
and could be very helpful for designing high magnetore- 
sistance graphene devices. 
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